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Abstract. The influence of the Pt cutoff on the production of bb events in Pythia 
has been studied. Two samples of events generated with P t > lGeV (default) and 
with P t > 10 GeV are investigated. In versions 5.75 and 6.115 of Pythia there are 
some critical differences between the two samples for what concerns contribution in bb 
production mechanisms. Version 6.152, taken as the CMS official one, has been tested 
also. 



1. Introduction 

CMS (Compact Muon Solenoid) is a general purpose experiment under construction at 
CERN. The description of the apparatus is given elsewhere CMS will be built on 
LHC (Large Hadron Collider) app collider with -^5=14 TeV and a maximum luminosity 
of 10 34 cm~~ 2 s^ 1 . One of the main physics topics of the Collaboration will be the 
study of the B hadron physics. In fact apart from the interest "per se" of this topic, 
b-jets are also the main background for several other physics channels ||, M. In this 
article, I present a study on bb production mechanisms at LHC energies using the Monte 
Carlo package Pythia || as event generator. 

At a hadron collider bb pairs are produced by the following mechanisms: 

(i) Pair production (gluon or quark fusion): gg — > bb, qq — > bb. For each event a bb 
pair is produced in the hard interaction. 

(ii) Flavour excitation: bq — ► bq, bg — > bg. A heavy flavour from the parton 
distribution of one beam particle is put on mass shell by scattering against a parton 
of the other beam. As b quark is not a valence flavour it must come from a branching 
g -> bb. 

(iii) Gluon splitting: g — > bb. It occurs in the initial- or final-state shower but no 
heavy flavours are produced in the hard scattering. 

A more detailed definition of the mechanisms can be found in ref ||. As these 
mechanisms give rise to different kinematic configurations of the final state is important 
to understand their contribution to the total cross section. 

The only way to generate (i), (ii) and (iii) in Pythia is using the steering card MSEL=1, 
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which generates QCD 2 — > 2 events. In order to have higher statistics, and save 
computation time the user can set kinematic cuts at generation level which limits the 
available portion of phase space J7), increasing the efficiency of selection cuts. One of 
this cuts is the P t cutoff. | 

In the following sections I will show differences between the use of the Pythia default 
cutoff (1 GeV) and a 10 GeV cutoff. Three versions of the generator have been 
considered: 5.75, 6.115 and 6.152 which has been chosen as event generator for the 
next CMS physics studies. 



2. bb production cross sections 



Two samples have been prepared with the three versions of Pythia to investigate the 
influence of the Pt cutoff on the production of bb events. Both of them have been 
generated using MSEL=1 and contain events with only one bb pair. The first sample 
(SAMPLE A) has been generated with the default P t cutoff (1 GeV). The second sample 
(SAMPLE B) has been generated with P t > 10 GeV. For both samples the bb production 
cross section for P t > 10 GeV has been computed using the following formula: 
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a 



(Pt > 10 GeV) 
bb 



N, 



bb 



&tot 



(1) 



tot 



where iV7 r '- e is the number of bb events with P t > 10 GeV, N tot is the total 

bb 1 — ' 

number of events generated and a tot is the total (QCD 2^2) cross section given by 
Pythia . The results of the computation for the three versions of the package are shown 
in Tab.|I] (the values of the cross section are given in /A>). 



cross section (P t > 10 GeV) (/ib) 


P 5.75 


P 6.115 


P 6.152 




A 


B 
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A 
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° pair production 


20 


20 


47 


51 


48 


52 


^flavour excitation 




110 




221 


211 


236 


°gluon splitting 


130 


130 


173 


193 


159 


173 


a bb 


150 


260 


220 


465 


418 


461 



Table 1. Contributions to the bb production cross section for P t > 10 GeV with 
Pythia 5.75, 6.115 and 6.152 (in fib). 



All double countings have been avoided requiring the number of b-quarks coming 
from the hard interaction to be 2 (pair production), 1 (flavour excitation) or (gluon 
splitting). As it appears clear from the table for version 5.75 and 6.115 the cross sec- 
tion is Pt dependent. The problem is related to flavour excitation which in this two 

|| Pt is defined as the transverse momentum of the outgoing partons in the center of mass system of 
the ingoing partons. 
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versions cannot be generated with default P t cutoff (sample A). The difference in the 
a b i (for sample A) between 5.75 and 6.115 is due to the different choice of the default 
parton distribution function. The interpretation of the observed phenomena, given by 
T. Sjostrand and E. Norrbin, is reported in ||. In version 6.152 the problem related to 
flavour excitation has been set but there is still a difference of about 10% between the 
cross sections of the two samples (for all versions statistical errors can be neglected). 
The difference in the value of cross sections is not very important, because the results 
are usually normalized to the total bb cross section of 500 /ib. On the other hand, the 
different contributions for the production mechanisms could be more dangerous, as they 
can lead to different kinematic distributions, and therefore affect the efficiencies for dif- 
ferent selections. 

Fig. [I] illustrates the difference in the bb production cross sections due to the additional 
contribution of the flavour excitation mechanism in sample B. 



cross section Vs pthat 




Figure 1. Cross section with the two different cutoff on P t . The solid line is with the 
default cutoff (1 GeV). The picture is obtained with Pythia 5.75. 



3. Kinematic distributions 

The main kinematic parameters which define the signature of bb event are the transverse 
momenta and pseudorapidities of the b quarks, and the angular distance A0 between 
their directions in the transverse plane. The first two parameters have similar 
distributions in both samples. The A<p distribution is shown in fig. |2] for the three 
different mechanisms. For what concerns gluon splitting, the distribution is slightly 
peaked at small A<p. The angle between the two 6-quarks produced by the gluon-fusion 
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Figure 2. Distribution of the angle between the two b quarks in the transverse plane: 
the upper one is gluon splitting, the middle is gluon fusion, the last is related to flavour 
excitation. The plots are in arbitrary units. 



mechanism has a peak at A</> ~ n, as expected, since in the process gg — > bb the b-quarks 
are produced back-to-back in the transverse plane. The last distribution corresponds 
to the flavour excitation production mechanism, for which the back-to-back topology 
is preferred. We can conclude that the total A<f> distributions of sample A and sample 
B are slightly different. Some care should be taken about this, as it could affect the 
estimated efficiency of selection cuts. 
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4. Conclusions 

Versions 5.75 and 6.115 show a big discrepancy in the total bb production cross section 
when using two different cutoffs on P t . The problem is due to flavour excitation which 
can be generated only in sample B but gluon splitting and pair production are stable 
and show the same cross section in the two samples. In both versions the use of a 
different P t . cutoff from the default one leads to slightly different A0 b ^ distributions 
caused by the presence of flavour excitation. 

In version 6.152 everything is set, a part a little inconsistency of about 10% on which I 
am working in order to understand the cause. 
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